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The cosmic acceleration is one of the most significant cosmological discoveries over the 
last century. Following the more accurate data a more dramatic result appears: the recent 
analysis of the observation data (especially from SNe la) indicate that the time varying dark 
energy gives a better fit than a cosmological constant, and in particular, the equation of 
state parameter w (defined as the ratio of pressure to energy density) crosses —1 at some 
low redshift region. This crossing behavior is a serious challenge to fundamental physics. In 
this article, we review a number of approaches which try to explain this remarkable crossing 
behavior. First we show the key observations which imply the crossing behavior. And then 
we concentrate on the theoretical progresses on the dark energy models which can realize 
the crossing —1 phenomenon. We discuss three kinds of dark energy models: 1. two-field 
models (quintom-like), 2. interacting models (dark energy interacts with dark matter), and 
3. the models in frame of modified gravity theory (concentrating on brane world). 
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I. THE UNIVERSE IS ACCELERATING 

Cosmology is an old and young branch of science. Every nation had his own creative idea 
about this subject. However, till the 1920s about the unique observation which had cosmological 
significance was a dark sky at night. On the other hand, we did not prepare a proper theoretical 
foundation till the construction of general relativity. Einstein's 1917 paper is the starting point of 
modern cosmology [l|. The next mile stone was the discovery of cosmic expansion, that was the 
recession of galaxies and the recession velocity was proportional to the distance to us. 

Except some rare cases, our researches are always based on the cosmological principle, which 
says that the universe is homogeneous and isotropic. In the early time, this is only a supposition 
to simplify the discussions. Now we have enough evidences that the universe is homogeneous and 
isotropic at the scale larger than 100 Mpc. The cosmological principle requires that the metric of 
the universe is FRW metric, 

ds 2 = -dt 2 + a 2 (t)(dr 2 + r 2 dn 2 2 ); (1) 
ds 2 = -dt 2 + a 2 (t)(dr 2 + sm(r) 2 dn 2 2 ); (2) 
ds 2 = -dt 2 + a 2 (t)(dr 2 + sinh(r) 2 (i^), (3) 

depending on the spatial curvature, which can be Euclidean, spherical or pseudo-spherical. Here, 
t is the cosmic time, a denotes the scale factor, r represents the comoving radial coordinate of the 
maximal symmetric 3-space, and dQ?, stands for a 2-sphere. Which geometry serves our space is 
decided by observations. FRW metric describes the kinetic evolution of the universe. To describe 
the dynamical evolution of the universe, that is, the function of a(t), we need the gravity theory 
which ascribes the space geometry to matter. The present standard gravity theory is general 
relativity. In 1922 and 1924, Friedmann found that there was no static cosmological solution in 

n 

general relativity, that is to say, the universe is either expanding or contracting [2]. To get a static 
universe, Einstein introduce the cosmological constant. However, even in the Einstein universe, 
where the contraction of the dust is exactly counteracted by the repulsion of the cosmological 
constant, the equilibrium is only tentative since it is a non-stationary equilibrium. Any small 
perturbation will cause it to contract or expand. Hence, in some sense we can say that general 
relativity predicts an expanding (or contracting) universe, which should be regarded as one of the 
most important prediction of relativity. 

In almost 70 years since the discovery of the cosmic expansion in 1929 [3], people generally 
believe that the universe is expanding but the velocity is slowing down. People try to understand 
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via observation that the universe will expand forever or become contracting at some stage. A 
striking result appeared in 1998, which demonstrated that the universe is accelerating rather than 
decelerating. Now we show how to conclude that our universe is accelerating. We introduce the 
standard general relativity, 

G F + Av = 8^GV, (4) 

where G^ v is Einstein tensor, G is (4-dimensional) Newton constant, ta V denotes the energy- 
momentum tensor, g^ u stands for the metric of a spacetime, fi, v run from to 3. Throughout 
this article, we take a convention that c = h = 1 without special notation. Define a new energy 
momentum T^ u , 

8ttGT^ = 8ttGV - Ag^. (5) 

T^y has included the contribution of the cosmological constant, whose effect can not be distin- 
guished with vacuum if we only consider gravity. 

The 00 component of Einstein equation dH) is called Friedmann equation, 

n k 8ttG , , 

"+;? = — * (6) 

where H = ^ is the Hubble parameter, an overdot stands for the derivative with respect to the 
cosmic time, k is the spatial curvature of the FRW metric, for ([T]), k = 0; for ([2]), k = 1; for 
([3]), k = —1, p = —Tq. Throughout this article, we take the signature (—,+,...,+). The spatial 
component of Einstein equation can be replaced by the continuity equation, which is much more 
convenient, 

p + 3H(p + p) = 0, (7) 

where p = Tj = T| = T| . Here we use a supposition that the source of the universe is in 
perfect fluid form. Using ([6]) and ([7j), we derive the condition for acceleration, 

a 4irG , „ . 

- = 5-(P + 3p). (8) 

We see that the universe is accelerating when p + 3p < 0. We know that the galaxies and dark 
matter inhabit in the universe long ago. They are dust matters with zero pressure. Thus, if 
the universe is accelerating, there must exist an exotic matter with negative pressure or we should 
modify general relativity. The cosmological constant is a far simple candidate for this exotic matter, 
or dubbed dark energy. For convenience we separate the contribution of the cosmological constant 
(dark energy) from other sectors in the energy momentum, the Friedmann equation becomes, 
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H 2 

2=n A0 + n m0 (i + z) 3 + n k0 (i + z) 2 , (9) 



#0 



where z is the redshift, the subscript denotes the present value of a quantity, 

°- = 8^G' (10) 



^kO ~ -77272- ( 12 ) 

H a 

The type la supernova is a most powerful tool to probe the expanding rate of the universe. 
In short, type la supernova is a supernova which just reaches the Chandrasekhar limit (1.4 solar 
mass) and then explodes. Hence they have the same local luminosity since they have roughly the 
same mass and the same exploding process. They are the standard candles in the unverse. We 
can get the distance of a type la supernova through its apparent magnitude. A sample of type la 
supernovae will generates a diagram of Hubble parameter versus distance, through which we get 
the information of the expanding velocity in the history of the universe. In 1998, two independent 
groups found that the universe is accelerating using the observation data of supernovae 
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that the data accumulate fairly quickly. The famous sample includes Gold04 [jj, Gold06 



4]. After 
fi|, SNLS 



ESSENCE [8j, Davis07 [9j, Union [101 ] . Constitution [ll| . Here, we show some results of one of 
the most recent sample, Union [lo| . which is plotted by x 2 statistics. 

From figHJ we see that: 1. the universe is almost spatially flat, that is the curvature term Q^o 
is very small. 2. the present universe is dominated by cosmological constant (dark energy), whose 
partition is approximately 70%, and the partition of dust is 30%. We introduce a dimensionless 
parameter, the deceleration parameter q, 

g = -p. (13) 

A negative deceleration parameter denotes acceleration. In a universe with dust and cosmological 
constant (which is called ACDM model), by definition 

q = - JJa, (14) 
whose present value q$ = —0.55. Hence the present universe is accelerating. 
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FIG. 1: The figure displays the counter constrained by SNe la (in blue), by BAO (in green), by CMB (in 
yellow), and the joint constraint by all the three kinds of observation (in black and white). £Ia is ^ao i n 
(flu]) . fl m is f2 m o in (fTTj) . This figure is borrowed from [id ]. 

The previous result depends on a special cosmological model, ACDM model in frame of general 
relativity. How about the conclusion if we only consider the kinetics of the universe? 

The simplest kinetic model is a sudden transition model, in which the deceleration parameter 
is a constant in some high redshift region and jumps to another constant at a critical redshift. The 
other simple choice is that the deceleration parameter is a linear function of z. We show the two 
kinetic models with the dynamic model ACDM in fig EJ It is clear that the universe accelerates in 
the present epoch in all the three models. A more rigorous analyze shows that the evidence for an 
accelerating universe is fairly strong (more than 5 a) [12!]. So we should investigate it seriously. 

ACDM is the most simple model for the acceleration, which is a concordance model of several 
observations. As we shown in fig[H the counters of CMB, BAO, and SNe la have cross section, which 
almost laps over result of the joint fittings. However, ACDM has its own theoretical problems. 
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FIG. 2: Kinetic universe vs dynamic universe. The fitting results of sudden transition model, linear expansion 
model, and ACDM model. The solid line is the best-fit of the sudden transition model (the deceleration 
parameter jumps at some rcdshift); the long-dashed line denotes the best-fit of linear expansion model 
(q = Qo + Qi z - ?i is a constant) the short-dashed line represents the best-fit of ACDM model. From [l^ . 

Furthermore, it is found that a dynamical dark energy model fits the observation data better. 
Especially, there are some evidences that the equation of state (EOS) of dark energy may cross 
— 1, which is a serious challenge to the foundation of theoretical physics. 

In the next section we shall study some problems of ACDM model and display that a dynamical 
dark energy model is favored by observations. We'll focus on the crossing behavior implied by the 
observation. In section III, we study 3 kinds of models with a crossing phantom divide dark energy. 
In section IV, we present the conclusion and more references of this topic. 



II. A DARK ENERGY WITH CROSSING -1 EOS IS SLIGHTLY FAVORED BY 



ACDM has two famous theoretical problems. 

The first is the finetune problem. The effect of the vacuum energy can not be distinguished 
from the cosmological constant in gravity theory. We can calculate the vacuum energy by a well- 
constructed theory, quantum field theory (QFT), which says that the vacuum energy should be 
larger than the observed value by 122 orders of magnitude, if QFT works well up to the Planck 
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The problems of ACDM 



7 



scale. In super symmetric (SUSY) theory, the vacuum energy of the Bosons exactly counteracts 
the vacuum energy of Fermions, such that we obtain a zero vacuum energy. However, SUSY must 
break at the electro-weak scale. At that scale, the vacuum energy is still large than the observed 
value by 60 orders of magnitude. So for getting a vacuum energy we observed, we should introduce 
a bare cosmological constant Ab are - The effective vacuum energy /Oeffect then becomes, 

Reflect = g - ^-bare ~\~ /^vacuum- (15) 

g^jjAbarc and /^vacuum have to almost counteract each other but do not exactly counteract each 
other, leaving a tiny tail which is smaller than the /Vacuum by 60 orders of magnitude. Which 
mechanism can realize such a miraculous counteraction? 

The second problem is coincidence problem, which says that the cosmological constant keeps a 
constant while the density of the dust evolves as (1 + z) 3 in the history of the universe, then why 
do they approximately equal each other at "our era" ? Different from the first problem, the second 
problem says the present ratio of dark energy and dark matter is sensitively depends on the initial 
conditions. Essentially, the coincidence problem is the problem of an unnatural initial condition. 
The densities of different species in the universe redshift with different rate in the evolution of the 
universe, so if their densities coincidence in our era, their density ratio must be a specific, tiny 
number in the early universe. It is also a finetune problem, but a finetune problem of the initial 
condition. 

Except the above theoretical problems, ACDM also suffers from observation problem, especially 
when faced to the fine structure of the universe, including galaxies, clusters and voids. Some 
specific observations differ from the predictions of ACDM (with standard partitions of dust and 
cosmological constant) at a level of 2a or higher. Six observations are summarized in [13|]: 1. 
scale velocity flows is much larger than the prediction of ACDM, 2. Type la Supernovae (Sne la) 
at High Redshift are brighter than what ACDM indicates, 3. the void seems more empty than 
what ACDM predicts, 4. the cluster haloes look denser than what ACDM says, 5. the density 
function of galaxy haloes is smooth, while ACDM indicates a cusp in the core, 6. there are too 
much disk galaxies than the prediction of ACDM. We do not fully understand the dynamics and 
galaxies and galaxy clusters, that is, the gravitational perturbation theory at the small scale. The 
agreement may approve when we advance our perturbation theory with a cosmological constant 
and the simulation methods at the small scale. However, in the cosmological scale, there are also 
some evidences that the dark energy is dynamical, including no. 2 of the previous 6 problems. 
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B. crossing — 1 

With data accumulation, observations which favor dynamical dark energy become more and 
accurate. Now we loose the condition that p = —p for the exotic matter (dark energy) which 
accelerates the universe. We go beyond the ACDM model. We permit that the EOS of dark 
energy is not exactly equal to —1, but still a constant. The fitting results by different samples 
of SNe la are displayed in fig El We see that although a cosmological constant is permitted, the 
dark energy whose EOS < — 1 is favored by SNe la. The essence whose EOS is less than —1 is 
called phantom, which can be realized by a scalar field with negative kinetic term. The action for 
phantom ip is 

S ph = J d A x^~g (Up m ~ W)) , (16) 

where yj—g is the determinate of the metric, the lowercase Greeks run from to 3, U{ip) denotes 
the potential of the phantom. In an FRW universe, the density and pressure of the phantom (|16p 
reduce to 

P = -^ 2 + U, (17) 

p=-^-U, (18) 
respectively. Now the EOS of the phantom w is given by, 



p _ -U 

P ~ -\i> 2 + u'' 



™ = F -= \\ ,1 , (19) 



¥ 2 

which is always less than —1 for a positive U. It seems that phantom is proper candidate for the 

rn 

dark energy whose EOS less than —1 [14]. It is very famous that a phantom field is unstable when 
quantized since the energy has no lower bound. It will transit to a lower and lower energy state. 
In this article we have no time to discuss this important topic for phantom dark energy. We would 
point out the basic idea for this issue, which requires the life time of the phantom is much longer 
than the age of the universe such that we still have no chance to observe the decay of the phantom, 
though it is fundamentally unstable. For references, see [3]. A completely regular quantum stress 
with w < — 1 is suggested in [if]]. 

If the dark energy really behaves as phantom at some low redshift region, it is an unusual 
discovery. But the dark energy may be more fantastic. In some model-independent fittings, the 
EOS of dark energy crosses —1, which is a really remarkable property and a serious challenge to 
our present theory of fundamental physics. 
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FIG. 3: The EOS of dark energy fitted by SNe la in a spatially flat universe, the contours display 68.3 %, 
95.4 % and 99.7% confidence level on w and fl m o (f2 m in the figure). The results of the Union set are shown 
as filled contours. The empty contours, from left to right, show the results of the Gold sample, Davis 07, 
and the Union without SCP nearby data. From [l(J. 
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FIG. 4: Panel (a): Uncorrelated band-power estimates of the EOS w(z) of dark energy by SNe la (Gold 
set Q). Vertical error bars show the 1 and 2-cr error bars (in blue and green, respectively). The horizontal 
error bars denote the data bins used in [3] . Panel (b) : The window functions for each bin from low redshift 
to high redshift. Panel (c): the likelihoods of w(z) in the bins from low redshift to high redshift. From [ljj]. 



Pioneer results of the crossing —1 of EOS of dark energy appeared in [171. Il8l|. Fig [4] illuminates 
that the EOS of dark energy may cross —1 in some low redshift. In fig [H the Gold04 data are 
applied, a uniform prior of 0.22 < $7 m o < 0.38 is assumed, and a spatially flat universe is the 
working frame. 

The perturbation of the dark energy will growth if its EOS is not exactly —1 in the evolution 
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FIG. 5: Constraints on the EOS of w(z) by WMAP3 [20] and Gold04 jjj. The light grey region denotes 2 
cr constraint, while the dark grey for la constraint. The left panel shows the constraint with dark energy 
perturbation, while the right displays the result without dark energy perturbation. From [19|. 



history of the universe. Hence to fit a model with dynamical dark energy with observation, the 
perturbation of the dark matter should be considered in principle. Such a study was presented 
in [ig| . in which a parametrization of the EOS of the dark energy with two constant Wq, w\ was 
applied, 



W = Wq + W\ 



l + Z 



(20) 



The result is shown in fig [5j We see from fig [5] that there is a mild tendency that the EOS of the 
dark energy cross —1. For a more general parametrization of EOS for dark energy, see 21]. 



With more and accurate data, the possibility of crossing —1 (phantom divide) seems a little 



more specific, see for example 22]. This crossing behavior is a significant challenge for theoretical 
physics. It was proved that the EOS of dark energy can not cross the phantom divide if 1. a 
dark energy component with an arbitrary scalar-field Lagrangian, which has a general dependence 
on the field itself and its first derivatives, 2. general relativity holds and 3. the spatially flat 
Friedmann universe 



n 



231 ] , for a more detailed proof, wee the appendix of [24] • Thus realizing such a 
crossing is not a trivial work. In the next section we investigate the theoretical progresses for this 
extraordinary phenomenon. 
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III. THREE ROADS TO CROSS THE PHANTOM DIVIDE 



23[]. Now that the 



To cross the phantom divide, we must break at least one of the conditions in 
dark energy behaves as quintessence at some stage , while evolves as phantom at the other stage, a 
natural suggestion is that we should consider a 2-field model, a quintessence and a phantom. The 
potential is carefully chosen such that the quintessence dominates the universe at some stage while 
the phantom dominates the universe at the other stage. It was invented a name for such 2-field 
model, "quintom" . There are also some varieties of quintom, such as hessence. We introduce 
these 2-field models in the first subsection. The next road is to consider an interacting model, in 
which the dark energy interacts with dark matter. The interaction can realize the crossing behavior 
which is difficult for independent dark energy. We shall study the interacting models in subsection 
B. The other possibility is that general relativity fails at the cosmological scale. The ordinary dark 
energy candidates, such as quintessence or phantom, can cross the phantom divide in a modified 
gravity theory. We investigate this approach in subsection C. 

A. 2-field model 



A typical 2-field model is the quintom model, which was proposed in 25], and was widely 



investigated later [26]. Generally , the action of a universe with quintom dark energy S is 

V 



S = I d 4 xV=g ( ^ + Atufl ) , (21) 



where R is the Ricci scalar, £ s tuff encloses all kinds of the stuff in the universe, for instance the 
dust matter, radiation, and quintom. At the late universe, the radiation can be negligible. So, 
often we only consider the dark energy, here quintom £ qu intomi and dust matter £<jm> 

^quintom = "fy^ + \d^d^ - W(<j>, %b) . (22) 

In (I22p the first term is the kinetic term of an ordinary scalar, the second term is the kinetic term 
of a phantom, and W((f), tp) is an arbitrary function of <j> and tp. In an FRW universe, the density 
and pressure of the quintom are 

P = \<p 2 -^ 2 + W, (23) 
V=\ft~ \tf ~ W. (24) 
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Hence, the EOS of the quintom w is 

in — .. . . . . . 

U 2 -U 2 + w' 



U 2 -U 2 -w , 



u> = — 1 requires 

(/» 2 = ft. (26) 

We see that in a quintom model, we do not require a static field (a field with zero kinetic term or 
a field at ground state) to get a cosmological constant. We only need that ip and (ft evolves in the 
same step, w < — 1 implies, 

ft - ft < 0, (27) 

if 

^ 2 - \ft + W > 0; (28) 

and 



12 -.2 



V> z > 0, (29) 



if 



^ 2 - \ft + W < 0. (30) 

([3lJj) yields an unnatural physical result, that is, the density of dark energy is negative. However, 
this is not as serious as the first glance, since we have little knowledge of the dark energy besides 
its effect of gravitation. Several evidences imply that we should go beyond the standard model of 
the particle physics when we describe dark energy. There are a few dark energy models permit 
density of the dark energy, or a component of it is negative (at the same time keep the total 



density positive), for example, see 



27 



281 ] . But, for a model with only two components, a dust and 
a quintom, it is difficult to set a negative density dark energy. In that case we need too much dust 
than we observed or a big curvature term. In the following text of this section, we only consider a 
dark energy with positive density. So w > — 1 implies, 

ft - > 0. (31) 

In summary, if the kinetic term of the quintessence dominates that of phantom, the quintom 
behaves as quintessence; else it behaves as phantom. We should select a proper potential to make 
quintessence and phantom dominate alternatively such that we can realize the crossing behavior. 



13 



A simple choice of the potential is that the quintessence and the phantom do not interact with 
each other, which requires, W((j>, t(i) = V{(p) + U(ip). The exponential potential is an important 
example which can be solved exactly in the quintessence model (a toy universe only composed by 
quintessence). In addition, we know that such exponential potentials of scalar fields occur naturally 
in some fundamental theories such as string/M theories. We introduce a model with such potentials 
in 29(, in which the potential V(4>, ip) is given by 

W{4>, VO = V{<j>) + Uty) = A^e~ x ^ + A^e-^, (32) 

where and are the amplitude of the potentials, k? = 8ttG, and are two constants. 
Since there is no direct couple between the quintessence and the phantom, the equations of motion 
of the quintessence and the phantom are two independent equations, 

+ 3H<j) + ^— = 0, (33) 



i> + 3H?P - — = 0. (34) 
dtp 

The continuity equation of the dust reads, 

Pdust + 3Hp dust = 0, (35) 

where Pdust denotes the density of the dust. The method of dynamical system has been widely 
used in cosmology. This method can offer a clear history of the cosmic evolution, especially the 
final states of the university. For applying this method, first we define the following dimensionless 
variables, 
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Label x$ y^ x<f, y^ z Stability 

K -x\ + x\ = 1 unstable 

P -^L ^(l + ^r) stable 

5 unstable 

F 1 unstable 

TO * g fT^J unstable 



TABLE I: The critical points, from 
the evolution equations (|3"3"l) - ([3"5j) become, 



-3a^ ( 1 + x\ - x\ - ^z 2 ) + X^-yl , (37) 



%'cf, 

y't = 3 y</> (y~ x l + x l + \ z2 - h-Y x <t> j . ( 38 ) 

„' _ o™ / 1 , „2 ,.2 1 A \ V^,5 



= -3aty ( 1 + x+ - xj, - -z j - X^—y^ , (39) 



y'i> = (y~ x l + x l + \ z2 - x ^^- x ^j > ( 40 ) 

a 7 = 3^ ( -zj + x\ + iz 2 - , (41) 

in which a prime denotes derivative with respect to In a. Generally, z in the above set will not be 
confused with redshift. The five equations in this system are not independent. They are constrained 
by Fridemann equation, 

R2 = Y Q^ + y -^ 2 + f7 + ^ust) , (42) 

which becomes 

x l + v\ - 4 + y% + z 2 = l. (43) 

with the dimensionless variables defined before. The critical points dwell at x'^ = y'^ = x'^ = y'^ = 
z' = 0. We present the result in table HI For detailed discussion of the critical points, see [29']. We 
would like to show a numerical example in which the EOS of the quintom crosses the phantom 
divide. Fig 6 illuminates that the EOS crosses —1. 

The previous quintom model includes two fields, which are completely independent and rather 
arbitrary. We can impose some symmetry in the quintom model. An interesting model with an 
internal symmetry between the two fields which work as dark energy is hessence [30]. Rather than 
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FIG. 6: The evolution of the effective equation of state of the phantom and normal scalar fields with W(<fi, a) 
for the case = 1. From 29 1. 



two uncorrelated fields, we consider one complex scalar field with internal symmetry between the 
real and the imaginary parts, 



<3? = </>i + icf) 2 , 



(44) 



with a Lagrangian density 

£hess = ~\ [(d^f + (d^*) 2 ] ~ Vfc **) = ~\ [ (d,0 2 ~ ] - V(0, (45) 



which is invariant under the transformation, 



<Pi - > <Pi cos a — i<p2 sm a, 
(>2 — ► —i(j>\ sin a + 4>2 cos a, 



(46) 
(47) 



if the potential is only a function of <J? 2 + (<J?*) 2 . For convenience, in (|45p we have introduced two 
new variables (£,0), 



£cosh6>, 4>2 = ^sinh^, 



which are denned by 



£ 2 = <j)\ - <t>l, coth#=-ji 



The equations of motion of £ and 9 are 



• 2 . dV 



Z + 3HZ + ZP + — = 0, 



(48) 



(49) 



(50) 
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fe + (2££ + 3Hf)9 = 0. 



(51) 



Clearly, £ and 8 couple to each other. The pressure and density of the hessence read, 



1 



(52) 



(53) 



respectively. The EOS of hessence, playing as dark energy, 



w 



(54) 



Qualitatively, hessence evolves as quintessence when £ 2 > £ 2 # 2 , while as phantom when £ 2 < £ 2 # 2 . 
The Lagrangian (145p does not include 6, hence the canonical momentum ir^ corresponding to the 
cyclic coordinate 9 are conserved quantities, 



(55) 



In an FRW universe, only exists. We define a conserved quantity Q which is proportional to 
o 



7T, 



3f2/ 



Q = a s i 



With this conserved quantity, the EOS becomes, 



(56) 



w 



C2_ Q 2 



Ic2 
2 



v(0 



(57) 



which is only a function of £. The Friedmann equations read as 



H 



8ttG 



(58) 



where pdust is the energy density of dust. The continuity equation of dust is (j35j) . The continuity 
equations of hessence are identical to the equations of motion (|50p and (|5ip . Then the system is 
closed and we present a numerical example in fig 7. Evidently, the EOS of hessence, playing the 
role of dark energy, crosses —1 at about a = 0.95 (z = 0.06). 

After the presentation of hessence model, several aspects of this model have been investigated, 
including to avoid the big rip 



lessence by recent observations 
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311, a ttr actor solutions for general hessence 



33] 



321 ], reconstruction of 



341 ]. and holographic hessence model 



dynamics of hessence in frame of loop quantum cosmology 
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B. interacting model 

Two-field model is a natural and obvious construction to realize the crossing —1 behavior of 
dark energy. However, there are two many parameters in the set-up, though we can impose some 
symmetries to reduce the parameters to a smaller region. One symmetry decrease one parameter, 
but we have little clue to impose the symmetries since we have no evidence in the ground labs. 

Interaction is a universal phenomenon in the physics world. An interaction term is helpful 
to cross the phantom divide. To illuminate this point, we first carefully analyze the previous 



observations which imply the crossing. Both the results of [18| and [191 ]. which are shown in figH] 
and fig respectively, are derived with a presupposition, that is, the dark energy evolves freely. 
In fact, what we observed is the effective EOS of the dark energy in the sense of gravity at the 
cosmological scale. When we suppose it evolves freely, we find that its EOS may cross the phantom 
divide. We can demonstrate for an essence with (local) EOS< —1, the cosmological effective EOS 
can cross —1 by aids of an interacting term. For the case with interaction, the continuity equation 
for dark energy becomes, 

p de + 3H(p de + p de ) = -r, (59) 

or 

p dc + 3H( Pdc + Pde + — ) = 0. (60) 

Here pd e is the density of dark energy, pde denotes the local pressure measured in the lab (if we 
can measure), V stands for the interaction term, and p e g = Pde + jjj is the effective pressure in the 
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cosmological sense. In a universe without expanding or contracting, H = 0, the interaction does 
no effect on the continuity equation, or energy conservation law, and thus does not yield surplus 
pressure [ti| . Two special cases are interesting: 1. is a constant, under which the interaction 
term contributes a constant pressure throughout the history of the universe. 2. SI f pd is a constant, 
under which the interaction term contributes a constant EOS in the history of the universe. In 
frame of a quintessence or phantom dark energy, the interaction term S jf Pd only shifts the EOS up 
or down by a constant distance in the w — z plane, without changing the profile of the curve of w. 
While the term ^ shifts the pressure, which can change the EOS significantly since the density 
Pde is a variable in the history of the universe. 

If the dark energy can couple to some stuff of the universe, the dark matter is the best candidate. 
Although non-minimal coupling between the dark energy and ordinary matter fluids is strongly 



37| , due to the unknown nature of the dark 



restricted by the experimental tests in the solar system 
matter as part of the background, it is possible to have non-gravitational interactions between the 
dark energy and the dark matter components, without conflict with the experimental data. The 
continuity equation for dust-like dark matter reads, 

Pdm + 3Hp dm = T. (61) 

Based on the previous discussion, we assume a most simple case 

T = H5p dm , (62) 

where 5 is constant 



38, 
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4dl | . This interaction term shifts a constant to the EOS of dark matter, 
that is, it is no longer evolving as (1 + z) 3 . We uniformly deal with quintessence and phantom, 
which are often labeled by X, with a constant EOS wx- So, the continuity equation of dark energy 
can be written as, 

Px + 3H(p x + w x px) = ~H5p dm . (63) 

Integrating (|61|) . we derive 

Pdm = Pdm0a~ 3+5 = Pdm0(l + ^) 3 ~ <5 - (64) 

Substituting to (f63|) . we reach 



px = Pxo(l + z) 3 ( 1+ ^) + Afano * [(1 + z)^ 1 ^ - (1 + zf- 5 } . (65) 

o + 3wx L J 
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FIG. 8: wx vs 6 under the condition gmr^j = 0. 



Only from the above equation, we can extract the effective EOS of the dark energy. To see this 
point, we make a short discussion. In a dynamical universe with interaction, the effective EOS of 
dark energy reads, 

PcS Pde + T/3H 



W d e 



1 + i dln Pde 



(66) 

is less than 0, 



Pde Pde 3<Zln(l + z) 

Clearly, if ^T^fe * s g rea ter than 0, dark energy evolves as quintessence; if d\n(i+z) 
it evolves as phantom; if ^^fe equals 0, it is just cosmological constant. In a more intuitionistic 
way, if pde decreases and then increases with respect to redshift (or time), or increases and then 
decreases, which implies that EOS of dark energy crosses phantom divide. So, some time we 
directly use the evolution of density of dark energy to describe the EOS of it. There is a more 
important motivation to use the density directly: the density is more closely related to observables, 
hence is more tightly constrained for the same number of redshift bins used 
The derivative of px with respect to (1+z) reads, 
dpx . . 



d{l + z) 



3(l+w x )pxo(l+z) 2+3wx +p dl 



m0 



5 + 3wx 



3(l + wx)(l + z) 



2+3w x 



(3 - S)(l + zf- s 
(67) 

If jfi+z) = a t some redshift z = z c , the effective EOS crosses —1. The result is illuminated by 
fig 8, in which we set z c = 0.3 as an example. This figure displays the corresponding wx when 
one fixes a 5, or vice versa if we require the EOS crosses -1 at z c = 0.3. This is an original figure 
plotted for this review article. 

Then the Friedmann equation reads, 

= "XCK-L TZ) - ■ f - | - 

'0 



M^ l + zf{1+wx) + \ + , W x 



5(1 + z 



Mi+w x ) 



+ 3wx(l + z 



,3-5 



(68) 
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FIG. 9: Constraints of (wx, S) by SNLS data at 68.3%, 95.4% and 99.7% confidence levels marginalized over 
fijfo with priors Q. X o = 0.72 ± 0.04 and S < 3. From 



where tlxo = k 2 pxo/(3Hq), and we have used ^dmo + ^xo = 1- Thus we need to constrain the three 
parameters 5, Oxo; w x- The constraint result by SNLS data is shown in fig 9. 5 = and wx = — 1 
are indicated by the horizontal and vertical dashed lines, which represent the non- interacting 
XCDM model and interacting ACDM model, respectively. 

From fig 8 and 9, we see that the observations leave enough space for the parameters (5, wx) 
to cross the phantom divide. 

In the previous interacting model, we consider a phenomenological interaction, which is put in 
"by hand". We should find a more sound physical foundation for the interactions. We will deduce 
an interaction term from the low energy limit of string/M theory in the scenario of the interacting 
Chaplygin gas model j^ . 

The Cha plyg in gas model was suggested as a candidate of a unified model of dark energy and 
dark matter [42j . The Chaplygin gas is characterized by an exotic equation of state 

p ch = -A/pch, (69) 

where A is a positive constant. The above equation of state leads to a density evolution in the 
form 



Pch = ]jA + ^, (70) 

where B is an integration constant. The attractive feature of the model is that it naturally unifies 
both dark energy and dark matter. The reason is that, from (170p . the Chaplygin gas behaves as 
dust-like matter at early stage and as a cosmological constant at later stage. 
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Though Chaplygin gas has such a nice property, it is a serious flaw when one studies the 
fluctuation growth in Chaplygin gas model. It is found that Chaplygin gas produces oscillations 
or exponential blowup of the matter power spectrum, which is inconsistent with observations 361 ] . 



So we turn to a model that the Chaplygin gas only plays the role of dark energy. To cross the 
phantom divide we consider a model in which the Chaplygin gas couples to dark mater. 

Although non-minimal coupling between the dark e nerg y and ordinary matter fluids is strongly 
restricted by the experimental tests in the solar system 37]], due to the unknown nature of the dark 



matter as part of the background, it is possible to have non-gravitational interactions between the 
dark energy and the dark matter components, without conflict with the experimental data. Thus, 
the observation constrain the only proper candidate to be coupled to Chaplygin gas is dark matter. 

We consider the original Chaplygin gas, whose pressure and energy density satisfy the relation, 
p c h = —A/p c h. By assuming the cosmological principle the continuity equations are written as 

p ch + SHjchPch = -r, (71) 

and 

Pdm + 3H^ dm p dm = T, (72) 

where the subscript dm denotes dark matter, and 7 is defined as 

7 = l + - = l+w, (73) 
P 

in which w is the parameter of the state of equation, and "fdm = 1 throughout the evolution of the 
universe, whereas 7 C ^ is a variable. 

r is the interaction term between Chaplygin gas and dark matter. Since there does not exist any 
microphysical hint on the possible nature of a coupling between dark matter and Chaplygin gas (as 
dark energy), the interaction terms between dark energy and dark matter are rather arbitrary in 



literatures 43J]. Here we try to present a possible origin from fundamental field theory for T. 

Whereas we are still lack of a complete formulation of unified theory of all interactions (including 
gravity, electroweak and strong), there at present is at least one hopeful candidate, string/M theory. 
However, the theory is far away from mature such that it is still not known in a way that would 
enable us to ask the questions about space-time in a general manner, say nothing of the properties 
of realistic particles. Instead, we have to either resort to the effective action approach which takes 
into account stringy phenomena in perturbation theory, or we could study some special classes of 
string solutions which can be formulated in the non-perturbative regime. But the latter approach 
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is available only for some special solutions, most notably the BPS states or nearly BPS states in 
the string spectrum: They seems to have no relation to our realistic Universe. Especially, there still 
does not exist a non-perturbative formulation of generic cosmological solutions in string theory. 
Hence nearly all the investigations of realistic string cosmologies have been carried out essentially 
in the effective action range. Note that the departure of string-theoretic solutions away from 
general relativity is induced by the presence of additional degrees of freedom which emerge in the 
massless string spectrum. These fields, including the scalar dilaton field, the torsion tensor field, 
and others, couple to each other and to gravity non-minimally, and can influence the dynamics 
significantly. Thus such an effective low energy string theory deserve research to solve the dark 
energy problem. There a special class of scalar-tensor theories of gravity is considered to avoid 
singularities in cosmologies in j^. The action is written below, 



' -R - \d^d^ + — L-L dm (£, di, q- x g^) 



(74) 



16ttG 2 ^ Y ' q((j)) 2 ' 
where G is the Newton gravitational constant, <j> is a scalar field, L dm denotes Lagrangian of matter 
, £ represents different matter degrees of matter fields, q guarantees the coupling strength between 
the matter fields and the dilaton. With action (I74p . the interaction term can be written as follow 



44] 



F = HpJ^L. (75) 
a ma 

Here we introduce new variable q(a)' = q(a)^ Wn ~ 1 ^ 2 , where o is the scale factor in standard FRW 
metric. By assuming 

q'{a) = g oe 3 / c (^m+Pe)/Pd m dina ! ^ 

where pdm and p^ are the densities of dark matter and the scalar field respectively, one arrive at 
the interaction term, 

T = 3Hc{p dni + P<t> ). (77) 

With this interaction form we study the equation set ([TTj) and ([72]) . Set s = — ln(l + z), 
r = 3Hc( Pch + Pdm ), u = {ml)- 1 ^ 2 )-^^, v = {m 2 )- 1 ^ 2 )-^^ A' = A(3H 2 )- 2 (3p 2 )- 2 , 
where c is a constant without dimension. Using these variables, (|7ip and (|72p reduce to 

^ = -3u + 3c(it + u), (78) 
^ = -3(v-A'/v)-3c{u + v). (79) 
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We note that the variable time does not appear in the dynamical system flTSj) and (|79j) because 
time has been completely replaced by redshift s = —ln(l + z). The critical points of dynamical 
system (|78|) and (l79l) are given by 



du dv , . 



The solution of the above equation is 



c 



u c = v c , (81) 

1 — c 

2 - (1 - c)A'. (82) 



We see the final state of the model contains both Chaplygin gas and dark matter of constant 
densities if the singularity is stationary. The final state satisfies perfect cosmological principle: 
the universe is homogeneous and isotropic in space, as well as constant in time. Physically T in 
(|72h plays the role of matter creation term C in the theory of steady state universe at the future 
time-like infinity. Recall that c is the coupling constant, may be positive or negative, corresponds 
the energy to transfer from Chaplygin gas to dark matter or reversely. A' must be a positive 
constant, which denotes the final energy density if c is fixed. Also we can derive an interesting and 
simple relation between the static energy density ratio 

(83) 



l + r s 
where 

r s = Jim 8*1. (84) 

z ~*~ l Pch 

To investigate the properties of the dynamical system in the neighbourhood of the singularities, 
impose a perturbation to the critical points, 

d{6u) 



ds 



-35u + 3c(5u + 5v), (85) 



^M. = -3(Sv + ^5v)- 3c(5u + Sv) . (86) 
ds vj. 

The eigen equation of the above linear dynamical system (8u, 5v) reads 

(A/3) 2 + (2 + — *— )A/3 + 2 - 2c 2 = 0, (87) 
1 — c 

whose discriminant is 

A = [(1 - c) 4 + (3/2 - c) 2 ]/(l - c) 2 > 0. (88) 
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FIG. 10: The plane v versus u. (a) left panel: We consider the evolution of the universe from redshift 
z = e 2 — 1. The initial condition is taken as u = 0, v — 400; u = 50, v = 350; u — 100, t> = 300; 
u = 120, w = 280 on the four orbits, from the left to the right, respectively. It is clear that there is a 
stationary node, which attracts most orbits in the first quadrant. At the same time the orbits around the 
neibourhood of the singularity is not shown clearly, (b) right panel: Orbit distributions around the node 



u c = v c c/(l - c), v c = y/(l -c)A'. From |27j 

Therefore both of the two roots of eigen equation (|87p are real, consequently centre and focus 
singularities can not appear. Furthermore only r s £ (0, oo), such that c G (0,1), makes physical 
sense. Under this condition it is easy to show that both the two roots of ([87)) are negative. Hence 
the two singularities are stationary. However it is only the property of the linearized system (|85|) 
and (|86l) . or the property of orbits of the neighbourhoods of the singularities, while global Poincare- 
Hopf theorem requires that the total index of the singularities equals the Euler number of the phase 
space for the non-linear system (|78|) and (|79|) . So there exists other singularity except for the two 
nodes. In fact it is a non-stationary saddle point at u = 0, v = with index —1. This singularity 
has been omitted in solving equations ([781) and ([79 j) . The total index of the three singularities is 
1, which equals the Euler number of the phase space of this plane dynamical system. Hence there 
is no other singularities in this system. From these discussions we conclude that the global outline 
of the orbits of this non-linear dynamical system (|78|> and (|79p is similar to the electric fluxlines of 
two negative point charges. Here we plot figs 10 and 11 to show the properties of evolution of the 
universe controlled by the dynamical system ([78)1 and (|79p . As an example we set c = 0.2, A' = 0.9 
in figs 10-12. 

Further, to compare with observation data we need the explicit forms of u(x) and v(x), especially 
v(x). We need the properties of in our model, which is contained in v(x), to compare with 
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FIG. 11: The plane v versus u. (a) left panel: To show the global properties of dynamical system (|78|) 
and (|79f we have to include some "unphysical " initial conditions, such as u = —100, v — —300, except for 
physical initial conditions which have been shown in figure ITTJ1 (b) right panel: Orbits distributions around 
the nodes. The two nodes u c = v c c/(l — c), v c = -\/(l — c)A' and u c — v c c/(l — c), v c = — -^/(l — c)A' 
keep reflection symmetry about the original point. Just as we have analyzed, we see that the orbits of this 
dynamical system are similar to the electric fluxlines of two negative point charges. From 27 1 



observations. Eliminate u(x) by using (178j) and (|79l) we derive 



__ + |1 + (1+ ^)/c]- + 3c„ + 3(l-c)(, + 



dv , ,, . . 
— + 3(v- A 1 v) 
ds 



/(3c) 



0, 



(89) 



which has no analytic solution. We show some numerical solutions in figure [T2l We find that for 
proper region of parameter spaces, the effective equation of state of Chaplygin gas crosses the 
phantom divide successfully. 

Up to now all of our results do not depend on Einstein field equation. They only depend on the 
most sound principle in physics, that is, the continuity principle, or the energy conservation law. 
Different gravity theories correspond to different constraints imposed on our previous discussions. 
Our improvements show how far we can reach without information of dynamical evolution of the 



universe. 



(I78p illuminates that the dark matter in this interacting model does not behaves as dust. Qual- 
itatively, the dark matter gets energy from dark energy for a positive c, and becomes soft, ie, its 
energy density decreases slower than (1 + z) 3 in an expanding universe. The parameter which 
carries the total effects of cosmic fluids is the deceleration parameter q. From now on we introduce 
the Friedmann equation of the standard general relativity. As a simple case we study the evolution 
of q in a spatially flat universe. So q reads 
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FIG. 12: v versus s. The evolution of v with different initial conditions u(— 2) = 0, v(— 2) = 400; u(— 2) = 
50, v(— 2) = 350; u(— 2) = 100, u(— 2) = 300 reside on the blue, red, and yellow curves, respectively. 
Obviously the energy density of Chaplygin gas rolls down and then climbs up in some low redshift region. 
So the Chaplygin gas dark energy can cross the phantom divide w = — 1 in a fitting where the dark energy 
is treated as an independent component to dark matter. From 27], this figure has been re-plotted. 



aa _ 1 fu + v - 3A'/v 2 
^ a 2 2 \ u + v 

and density of Chaplygin gas u and density of dark matter v should satisfy 

u(0) +v(0) = 1. (91) 

And then Friedmann equation ensures the spatial flatness in the whole history of the universe. 
Before analyzing the evolution of q with redshift, we first study its asymptotic behaviors. When 
z — * 00, q must go to 1/2 because both Chaplygin gas and dark matter behave like dust , while 
when z — > — 1 q is determined by 

1 (u c + v c -3A' Ivl 
hm q = - — — £■ 

z-*-l 2 \ Uc + Vc 

One can finds the parameters c = 0.2, A' = 0.9 are difficult to content the previous constraint 
Friedmann constraint (|9ip . Here we carefully choose a new set of parameter which satisfies Fried- 
mann constraint (|9ip . say, A' = 0.4, c = 0.06. Therefore we obtain 

lim q = -1.95, (93) 

Z^r — 1 

by using ([HI]) and (f82j) . Then we plot figure [131 to clearly display the evolution of q. One can 
check u(0) = 0.25, v(0) = 0.75; u(0) = 0.28, u(0) = 0.72; u(0) = 0.3, u(0) = 0.7, respectively on 
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FIG. 13: q versus s. The evolution of q with different initial conditions u(— 2) = 0, v(— 2) = 273; u(— 2) = 
15, v(— 2) = 250; it(— 2) = 25, w(— 2) = 233, reside on the blue, red, and yellow curves, respectively. 
Evidently the deceleration parameter q of Chaplygin gas rolls down and crosses q = in some low redshift 
region. The transition from deceleration phase to acceleration phase occurs at z — 0.18; z = 0.21; z = 0.23 
to the curves u{-2) = 0, u(-2) = 273; u{-2) = 15, v(-2) = 250; u(-2) = 25, v(-2) = 233, respectively. 
One finds — g « 0.5 ~ 0.6 at z = 0, which is well consistent with observations. From [^l, this figure has 
been re-plotted. 

the curves v(— 2) = 273; u(— 2) = 250; u(— 2) = 233. One may find an interesting property of the 
deceleration parameter displayed in fig|13t the bigger the proportion of the dark energy, the smaller 
the absolute value of the deceleration parameter. The reason roots in the extraordinary state of 
Chaplygin gas (|69p , in which the pressure p c h is inversely proportional to the energy density . 

Also we note that maybe an FRW universe with non-zero spatial curvature fits deceleration 
parameter better than spatially flat FRW universe. This point deservers to research further. 

After the presentation of the original interacting Chaplygin gas model, there are several gen- 
eralizations. For details of these generalizations, see [3]. 




The judgement that there exists an exotic component with negative pressure, or dark energy, 
which accelerates the universe, is derived in frame of general relativity. The validity of general 
relativity has been well tested from the scale of millimeter to the scale of the solar system. Beyond 
this scale, the evidences are not so sound. So we should not be surprised if general relativity fails at 
the scale of the Hubble radius. Surely, any new gravity theory must reduce to general relativity at 




C. model in frame of modified gravity 
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the scale between millimeter to the solar system. In frame of the new gravity theories, the cosmic 
acceleration may be a natural result even we only have dust in the universe. 

There are various suggestions on how to modify general relativity. In this brief review we 
concentrates on the brane world theory. Inspired by the developments of string/M theory, the idea 
that our universe is a 3-brane embedded in a higher dimensional spacetime has received a great 
deal of attention in recent years. In this brane world scenario, the standard model particles are 
confined on the 3-brane, while the gravitation can propagate in the whole space. In this picture, 
the gravity field equation gets modified at the left hand side (LHS) in while the dark energy is a 
stuff put at the right hand side (RHS) in ((4]). In the modified gravity model, the surplus geometric 
terms respective to the Einstein tensor play the role of the dark energy in general relativity. 

We consider a 3-brane imbedded in a 5-dimensional bulk. The action includes the action of the 
bulk and the action of the brane, 

S = 5bulk + SWane- (^4) 

Here 

Sbuik = / d 5 ^v / -55^buik, (95) 
Jm 

where X = (t, Z . X . X . X 3 ) is the bulk coordinate, a; 1 ,^ 2 ,^ 3 are the coordinates of the maximally 
symmetric space. M denotes the bulk manifold. The bulk Lagrangian can be 

A>uik = ^[R 5 + aF(R 5 )]+£ m + A 5 , (96) 

where g$, K5, R$, C m , denote the bulk manifold, the determinant of the bulk metric, the 5- 
dimensional Newton constant, the 5-dimensional Ricci scalar, and the bulk matter Lagrangian, 
respectively. F(R$) denotes the higher order term of scalar curvature R§, the Ricci curvature 
i?5AB, the Riemann curvature -Rsabcd- 

There are too much possibilities and rather arbitrary to choose the higher order terms. Generally 
the resulting equations of motion of such a term give more than second derivatives of metric and 
the resulting theory is plagued by ghosts. However there exists a combination of quadratic terms, 
called Gauss-Bonnet term, which generates equation of motion without the terms more than second 
derivatives of metric and the theory is free of ghosts [46]. Another important property of Gauss- 
Bonnet term is that, just like Hilbert Lagrangian is a pure divergence in 2 dimensions and Einstein 
tensor identifies zero in 1 and 2 dimensions, we have that in 4 or less dimension the Gauss-Bonnet 
Lagrangian is a pure divergence. We see the dilemma of quadratic term in 4 dimensional theory: 
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if we include it with non pure divergence we shall confront ghosts; if we want to remove ghosts 
we get a pure divergence term. So only in theories in more than 4 dimensional Gauss-Bonnet 
combination provides physical effects. Moreover the Gauss-Bonnet term also appears in both low 
energy effective action of Bosonic string theory [47J and low energy effective action of Bosonic 
modes of heterotic and type II super string theory 48J. An investigation into the effects of a 
Gauss-Bonnet term in the 5 dimensional bulk of brane world models is therefore well motivated. 
The Gauss-Bonnet term in 5 dimension reads, 



F(i? 5 ) = Rl ~ 4#5AB ^ B + i?5ABCD ^ BC ° • (97) 

The action of the brane can be written as, 

Sbrane = / d^Xy/^g (k^ 2 K + L b rane) , (98) 
JM 

where M indicates the brane manifold, g denotes the determinant of the brane metric, Lbrane stands 
for the Lagrangian confined to the brane, and K marks the trace of the second fundamental form 
of the brane. x = (r, X j X j X ) is the brane coordinate. Note that r is not identified with t if the 
the brane is not fixed at a position in the extra dimension z =constant. We will investigate the 
cosmology of a moving brane along the extra dimension z in the bulk, and such that r is different 
from t. 

We set the Lagrangian confined to the brane as follows, 

ibrane = R ~ X + L m , (99) 

where A is the brane tension and L m denotes the ordinary matter, such as dust and radiation, 
located at the brane. R denotes the 4 dimensional scalar curvature term on the brane, which is 
an important one except a Gauss-Bonnet term in the bulk. This induced gravity correction arises 
because the localized matter fields on the brane, which couple to bulk gravitons, can generate via 
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quantum loops a localized four-dimensional world-volume kinetic term for gravitons 

Assuming there is a mirror symmetry in the bulk, we have the Friedmann equation on the brane 
[50! . see also 



51] 



' c 



1 + ( H 2 + — 



/.• r x ' x ' «■■-'■* N 2 



a 2 2 



where 




<'-±±^l> + »lT + %£)- (101) 
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r c = 4li 2 . (102) 

iere M is a constant, standing for the mass of bulk black hole. For various limits of (llOOp . see 
52]. 

For convenience, we introduce the following new variables and parameters, 

u = 3ff2 (p + A) = si m0 (i + z) 6 + n x , 

m = -aH , 
1 

n 



n r c 



l TTTT 2 1 ( r 4oAs 8aM( 5 )G\ 



2 " 3m \ V 6 a 

1 



w 1 + X /1 + ^A 5 +mO M0 (l + z) 4 ), (103) 

and we have assumed that there is only pressureless dust in the universe. As before, we have used 
the following notations 



k _8nGp m0 _8vrGA 3A 5 3Mk| 

^ A -^"^f' ^-^f' nM °-8^pf (10l) 



With these new variables and parameters, (|100|) can be rewritten as 

4n(x - 2y)[l + m(x + y)} 2 = (x - u) 2 . (105) 

This is a cubic equation of the variable x. According to algebraic theory it has 3 roots. One 
can explicitly write down three roots. But they are too lengthy and complicated to present here. 
Instead we only express those three roots formally in the order given in Mathematica 

x\ = x\(y, u\m, n), 
X2 = X2(y,u\m,n), 

X3 = x 3 {y,u\m,n), (106) 

where y and u are two variables, m and n stand for two parameters. The root on x of the equation 
(|105p gives us the modified Friedmann equation on the Gauss-Bonnet brane world with induced 
gravity. From the solutions given in (|106p . this model seems to have three branches. In addition, 
note that all parameters introduced in (|103j) and (|104p are not independent of each other. According 
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to the Friedmann equation (|106p . when all variables are taken current values, for example, z = 0, 
the Friedmann equation will give us a constraint on those parameters, 

1 = /(Ofco, Cl m0 , ^MOi ^A 5 , ^A, "1, n). (107) 

To compare with observation, we introduce the concept "equivalent dark energy" or "virtual 
dark energy" in the modified gravity models, since almost all the properties of dark energy are 
deduced in the frame of general relativity with a dark energy. 

The Friedmann equation in the four dimensional general relativity can be written as 



, k 8ttG . 

H 2 + - = —(p + Pde ), (108) 

Ob O 



where the first term of RHS of the above equation represents the dust matter and the second term 
stands for the dark energy. Generally speaking the Bianchi identity requires, 

^ + 3H( Pde + Pde ) = 0, (109) 

we can then express the equation of state for the dark energy as 

p de ldlnp de 
w d e = — = -!--— - . (110) 

Pde 3 alna 

Note that we can rewrite the Friedmann equation (|106[) in the form of (|108|) as 

xHl = -3-P+ ( H$x(y,u\m,n)-—p \ =-g-(f + Q), (HI) 

where p is the energy density of dust matter on the brane and the term 

3H 2 

Q = ^^x(y,u\m,n) - p (112) 

corresponds to pd e in (|108l) . 

In fig[T3]we show the equation of state for the virtual dark energy when we take m = 1.036 
and n = 0.04917. In this case, from the constraint equation (|128|) . one has CI mo = 2.08. From the 
figure we see that w e jf < — 1 at z = and 

dw ef f 



dz 



< 0. 

2=0 



Therefore the equation of state for the virtual dark energy can indeed cross the phantom divide 
w = — 1 near z ~ 0. 

Fig [H] illuminates that the behavior of the virtual dark energy seems rather strange. However 
we should remember that it is only virtual dark energy, not actual stuff. The whole evolution of the 
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W e ff 




m=l . 036 
n=0 .04917 



1.5 



2.5 



1 + z 



FIG. 14: The equation of state w e ff with respect to the red shift 1 + z, with f2 m o = 0.28 and Ha — 2.08. 
From [52| . 




1 + z 



FIG. 15: H 2 /H§ versus 1 + z, with n rnQ = 0.28 and il A = 2.08. From 



3- 



universe is described by the Hubble parameter. We plot the Hubble parameter H corresponding 
to figQIlin figCEl 

Fig [IS displays that the universe will eventually becomes a de Sitter one. For more figures with 
different parameters, see [52]. The constraint of this brane model with induced scalar term on the 



53] 



brane and Gauss-Bonnet term in the bulk has been investigated in 

The above is an example of "pure geometric" dark energy. We can also consider some mixed 
dark energy model, ie, the cosmic acceleration is driven by an exotic matter and some geometric 
effect in part. Why such an apparently complicated suggestion? There are many interesting 
models are proposed to explain the cosmic acceleration, including dark energy and modified gravity 
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models. However, several influential and hopeful models, such as quintessence and DGP model, 
fundamentally can not account for the crossing —1 behavior of dark energy. By contrast, some 
hybrid model of the dark energy and modified gravity may realize such a crossing. As an example 



54]. 



we study the quintessence and phantom in frame of DGP 

Our starting point is still action (|94p . In a DGP model with a scalar, (|96H becomes a pure 
Einstein-Hilbert action, 



£buik = TT - J^ 5 ' (113) 
Zk 5 



and we add a scalar term in (I99p . 



-^branc — ^Q n Q^ ^ m -^scalar- (H4) 

Here the scalar term can be ordinary scalar (quintessence) or phantom (scalar with negative kinetic 
term). The Lagrangian of a quintessence reads, 

L^—d^d^-VW, (115) 



and for phantom, 



^ = \d^d^ - U(i(>). (116) 



In an FRW universe we have 



P<t> = \<i> 2 + V(<t>), (117) 

P<t> = \j> 2 -V{4>). (118) 

The exponential potential is an important example which can be solved exactly in the standard 
model. Also it has been shown that the inflation driven by a scalar with exponential potential 
can exit naturally in the warped DGP model [551 ] . It is therefore quite interesting to investigate a 
scalar with such a potential in late time universe on a DGP brane. Here we set 

V = V e~ Xl i. (119) 

Here Ai is a constant and V$ denotes the initial value of the potential. 
The Friedmann equation (jlOOp becomes 



tjO k 1 

H + ~~2 = o 2 



p + p + 0p o (l + 

Po 



(120) 
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where 



p = (121) 

c 

Similar to the previous case, we derive the virtual dark energy by comparing (I120|) and (I108p . 



Pde = Pa + Po + 9p 



p + p o + 0po(l + ^) 1/2 



(122) 



From (lllOp . we calculate the derivation of effective density of dark energy with respective to 
ln(l + z) for a ordinary scalar, 

dpde or 12 | /a/-, , <A 2 + 2^ + 2p dm x m - 2 , 10 „> 

3[0 + (9(1 H ) ' (0 + pdm)\- (123) 



dln(l + z) po 

If = 1, both terms of RHS are positive, hence it never goes to zero at finite time. But if 6 = — 1, 
the two terms of RHS carry opposite sign, therefore it is possible that the EOS of dark energy 
crosses phantom divide. In a scalar-driven DGP, we only consider the case of 9 = — 1. 
For convenience, we define some dimensionless variables, 

» ~ vfe (124) 



» * ^ (125) 



A V Pm f-\na\ 
m = ^H> (126) 



y4 = Tl^- (127) 



The Friedmann equation (|120p becomes 



y? + y 2 2 + y 3 2 + y 2 - vl ( i + 2 ^ + | + yi ) ^ = 1. (128) 



The stagnation point, that is, dpd e / 'dln(l + z) = dwells at 



V4 (n,V3 



2 + ^f =2, (129) 



+ V ' vl, 

which can be derived from fj 123|) and (|128|) . 

One concludes from the above equation that a smaller r c , a smaller £l m o (Recall that it is 
defined as the present value of the energy density of dust matter over the critical density), or a 
larger (which is defined as the present value of the kinetic energy density of the scalar over 
the critical density) is helpful to shift the stagnation point to lower redshift region. We show a 
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FIG. 16: For this figure, fl^i — 0.01, Q rc = 0.01, Ai = 0.5. (a) The left panel: (3 and 7 as functions of s, 
in which (3 resides on the solid line, while 7 dwells at the dotted line. The EOS of dark energy crosses — 1 
at about s = —0.22, or z — 0.25. (b) The right panel: the corresponding deceleration parameter, which 



crosses at about s 



-0.40, or z = 0.49. From 
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concrete numerical example of this crossing behaviours in fig [16] For convenience we introduce the 
dimensionless density and rate of change with respect to redshift of dark energy as below, 



Pde 



Pc b z 



(130) 



where p c denotes the present critical density of the universe, and 

1 v\ dpde 



1 



p c VL Tc ds 



(1 + 2 y i +y l +y * )-^(2yl + yl) - 2y\ 



vl 



(131) 



A significant parameters from the viewpoint of observations is the deceleration parameter q, 
which carries the total effects of cosmic fluids. We plot q in these figures for corresponding density 
curve of dark energy. In the fig [16] we set £l m = 0.3. £l rc is defined as the present value of the 
energy density of po over the critical density Q rc = po/p c - 

Now, we turn to the evolution of a universe with a phantom (I116P in DGP. In an FRW universe 
the density and pressure of a phantom can be written as (fT7|) . (fl~8]) . 

To compare with the results of the ordinary scalar, here we set a same potential as before, 



U = U n e 



(132) 



The ratio of change of density of virtual dark energy with respective to ln(l + z) becomes, 



dln(l + z) 



3[-i; 2 + 9(l + 



-ip 2 + 2U + 2p dm x 



Po 



)-^\-^ + Pdm )]. 



(133) 
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To study the behaviour of the EOS of dark energy, we first take a look at the signs of the terms 
of RHS of the above equation. (— ip 2 + pdm) represents the total energy density of the cosmic 
fluids, which should be positive. The term (1 + -V 1 ' +zu+2p dm ^ — 1/2 gho^d a i so b e positive. Hence if 



1, both terms of RHS are negative: it never goes to zero at finite time. Contrarily, if 9 = 1, 



the two terms of RHS carry opposite sign: the EOS of dark energy is able to cross phantom divide. 
In the following of the present subsection we consider the branch of 6 = 1. 
Now the Friedmann constraint becomes 

- Vl + Vl + Vl + Vl + Vl (l + 2 - y * + ^ + y2 ^ = 1. (134) 

Again, one will see that in reasonable regions of parameters, the EOS of dark energy crosses —1, 
but from below —1 to above —1. 

The stagnation point of pd e inhabits at 

which can be derived from (|133p and (|134|) . One concludes from the above equation that a smaller 
r c , a smaller fi m , or a larger £1^ is helpful to shift the stagnation point to lower redshift region, 
which is the same as the case of an ordinary scalar. Then we show a concrete numerical example 
of the crossing behaviour of this case in fig [TTJ The dimensionless density and rate of change with 
respect to redshift of dark energy become, 

Pde SI 



and 



P Pc V 



7 



-vl + vl + vl + yf (i + 2 01+^ + ^ 1/2 



1 + 2 y l±4 +yl r 1/2 (-2yj + Vl) + 2y\ 

vi 



(136) 



(137) 



Similarly, the deceleration parameter is plotted in the figure [T71 for corresponding density curve of 
dark energy. In this figures we also set il m = 0.3. 

Fig 1171 explicitly illuminates that the EOS of virtual dark energy crosses —1, as expected. At 
the same time the deceleration parameter is consistent with observations. 



IV. SUMMARY 



The recent observations imply that the EOS of dark energy may cross —1. This is a remarkable 
phenomenon and attracts much theoretical attention. 
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FIG. 17: For this figure, Q, ki = 0.01, tt rc = 0.01, A = 0.01. (a) The left panel: (3 and 7 as functions of s, 
in which (3 resides on the solid line, while 7 dwells at the dotted line. The EOS of dark energy crosses — 1 
at about s = —1.25, or z = 1.49. (b) The right ganel: The corresponding deceleration parameter, which 



crosses at about s 



-0.50, or z = 0.65. From 
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We review three typical models for the crossing behavior. They are two-field model, interacting 
model, and modified gravity model. 

There are several other interesting suggestions in or beyond the three categories mentioned 
above. We try to list them here for the future researches. We are apologized for this incomplete 
reference list on this topic. 

Almost in all dark energy models the dark energy is suggested as scalar. However, 3 orthogonal 
vectors can also play this role. For the interacting vector dark energy and phantom divide crossing, 
sec 



561 ] . For the suggestion of crossing the phantom divide with a spinor, see [57]]. Multiple k- 
essence sources are helpful to fulfil the condition for phantom divide crossing [o^ ] . Phantom divide 
crossing can be realized by non-minimal coupling and Lorentz invariance violation 5j|. An exact 



60] 



s with crossing —1, see 



62]. 



fill]. The 



solution of a two-field model for this crossing has been found in 

For previous interacting X (quintessence or phantom) modei 
cosmology of interacting X in loop gravity has been studied in 

Interacting holographic dark energy is a possible mechanism for the phantom divide crossing 
631 ] . And the thermodynamics of interacting holographic dark energy with phantom divide crossing 



is investigated in 



641 ] . An explicit model of F(R) gravity in which the dark energy crosses the 



651 ] . The phantom-like effects in a DGP-inspired F(R, 



phantom divide is reconstructed in 
gravity model is investigated in 
space 67j, a new braneworld in the sourced-Taub background is proposed 681 ] . while the previous 



661 ] . Based on the recent progress in studies of source of Taub 



brane world models are imbedded in AdS (RS) or Minkowski (DGP). In this model the EOS for 
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the virtual dark energy of a dust brane in the source region can cross the phantom divide. For 
other suggestions in brane world model, see 69]. 

Similar to the coincidence problem of dark energy, we can ask why the EOS crosses —1 recently? 



This problem is studied in 



ZC]. 



On the observational side, the present data only mildly favor the crossing behavior. We need 
more data to confirm or exclude it. 

Theoretically, we should find more natural model which has less parameters. We must go beyond 
the standard model of particle physics. The problem cosmic acceleration is a pivotal problem to 
access new physics. To study the problem of crossing —1 EOS will impel the investigation to the 
new Laws of nature. 
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